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D I Abstract. We rederive in a simplified version the Lehmann-Sommers eigenvalue 

distribution for the Gaussian ensemble of asymmetric real matrices, invariant under 
real orthogonal transformations, as a basis for a detailed derivation of a Pfaffian 
generating functional for n-point densities. This produces a simple free-fermion 
• diagram expansion for the correlations leading to quaternion determinants in each 

order n. All will explicitly be given with the help of a very simple symplectic kernel 
for even dimension N. The kernel is valid both for complex and real eigenvalues 
and describes a deep connection between both. A slight modification by an artificial 
additional Grassmannian solves also the more complicated odd- A case. As illustration 
g ■ we present some numerical results in the space C™ of complex eigenvalue n-tuples. 
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O ! 1. Introduction 
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In a recent short communication [1] a simple derivation of the n-point eigenvalue cor- 
relations for the real Ginibre ensemble was announced, which will be presented here 
in more detail and extended to most general cases. Ginibre [2], when he proposed his 
three types of Gaussian non-Hermitian matrix ensembles with complex, quaternion real 
and real entries respectively, was not able to solve the correlations for the real ensemble 
invariant under orthogonal transformations. It took quite a long time until Lehmann 
and Sommers [3] derived the joint probability density of eigenvalues, which is somewhat 
difficult to understand since the eigenvalues can be real or pairwise complex conjugate. 
Below we will present a simplified version of the derivation since we need it for obtain- 
ing a generating functional for the correlations, which we present as symmetric n-point 
densities R n (zi, z 2 , ■ ■ ■ , z n ) in the space of complex eigenvalue n-tuples (z\, z 2 , ■ ■ ■ , z n ). 
These include, as singular contributions, correlations with real eigenvalues and with 
complex conjugate pairs. 



X 
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Edelman [I] rederived the joint probability density of eigenvalues some years later and 
determined the 1-point density of complex eigenvalues. In a somewhat earlier paper 
Edelman, Kostlan and Shub [5] calculated the density of real eigenvalues and determined 
their probability. We will use both results here to calculate a very simple skew-symmetric 
kernel K.n(z\, z^) which governs all correlation functions. This kernel follows by a simple 
argument from Edelman's complex density, it helps us to invert a large antisymmetric 
matrix, the result of which turns out to be simply tridiagonal. This is of course possible 
only for even dimension N. In this paper we will also show how to generalize the result 
for odd N. In more recent work Akemann and Kanzieper calculated the smooth com- 
plex correlations as Pfaffians [6] , Sinclair derived a generating Pfaffian functional [7j and 
Forrester and Nagao [8] were able to determine the real-real and complex-complex corre- 
lations as Pfaffians with the help of skew-orthogonal polynomials. Without mentioning 
the paper [1] Borodin and Sinclair [9J generalized this to more general and also to crossed 
correlations. The general solution for odd dimension N is not contained in these papers. 

Here we will use Grassmannians to derive a simple free fermion zero-dimensional field 
theory, which presents the n-point density as a diagram expansion, that turns out to 
be a quaternion determinant of a 2n x 2n selfdual matrix equivalent to a Pfaffian of 
a certain skew-symmetric matrix. The fermionic Wick theorem helps us to work with 
the complicated algebra of the Pfaffians. The odd- dimensional case can be managed 
by introducing an additional artificial Grassmannian and again the resulting large an- 
tisymmetric matrix, which we need to invert for perturbation expansion, turns out to 
have a simple tridiagonal inverse. Finally we will present the correlation functions in a 
form analytic in the dimension N valid for even and odd N. As illustration we show 
numerical simulations for 1-, 2- and 3- point densities and numerical presentations of 
analytical results for 1- and 2- point correlation functions. 

Let us mention that the real asymmetric Gaussian ensemble has many applications in 
physics and social sciences, such as biological webs [10J, neural networks [TT] directed 
quantum chaos [12], financial markets [13] and quantum information theory [II] . Our 
paper is concerned with the correlations of the eigenvalues z% = A< of such real matrices 
Jij (1 ^ i,j ^ N) which fulfill the characteristic equation det(Jy — z) = and therefore 
are real or pairwise complex conjugate. 

2. Joint density of eigenvalues 

We recall the derivation of the joint density of eigenvalues for the real Ginibre ensemble, 
which we simplify in the following. The normalized measure for the A-dimensional real 
asymmetric random matrix is given by: 




(1) 
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For the joint density of eigenvalues Aj we first consider the case of dimension N = 2, 
because this already shows the essential features. 



2.1. Dimension N = 2 

'A B 



The matrix J = = OJO' 1 may be brought by an orthogonal transformation O to 

the form: 



J 



cos <p 
— sin i 



sin ( 

cos< 



Ai 6 
-6 A 2 



cos « 
sin < 



- sin - 

cos 6 



(A x + A 2 )/2 5 

-5 (Ai+A 2 )/2 



+ 



Ai — A 2 / cos 20 — sin 20 
— sin 20 — cos 20 



• (2) 



Here the range of parameters is Ai > A 2 , ^ ^ ir, — oo < 5 < oo. Ai,A 2 are the 
eigenvalues of the symmetric part of J, is the angle of rotation which diagonalizes the 
symmetric part of J and 5 determines the skew-symmetric part of J. Note that (° 50 ) 
commutes with f cos< t sm ^ ,) and thus cannot be brought to diagonal form by such a 

\— Bin 9 cos<p/ o o 

transformation. With the skew-symmetric matrix 0~ 1 dO = (_° d ^ ) we obtain 

dJ = 0(0~ 1 dOJ ' -JO~ l dO + dJ)0~ 1 . (3) 
For the Jacobian we may drop the rotation of the increment^ and from 



O^dOJ-JO^dO+dJ 



d\x d5 + o?0(A 2 — \i] 

-d8 + c?0(A 2 — Ai) d\ 2 



we find 



d{A,B,C,D) 
d((f>,5, A 1; A 2 ) 



A 2 - Ai A 2 - Ai 

1 -10 

10 

1 



2(Ax-A 2 ). 



Thus the measure in the new parametrization is given by 

1 



djj,(J) 



(2rr) s 



■ d6 dXt d\ 2 2(Ai - A 2 )e-( A ? +A ^ 252 )/ 2 



which is positive for Ai > A 2 , ^ ^ tt, — oo < 5 < oo. 
A check of normalization shows that J dfi(J) = 1. 



(4) 



(5) 



(6) 



Let us now go to the eigenvalues A± of J which are obtained from 

A 2 - (Ai + A 2 )A + A X A 2 + 5 2 = 0. 



Ai — A 5 

-5 A 2 - A 



Therefore 



(7) 



(8) 
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A, + A_ //A. - A \ 2 ro 



2 V V 2 

with Ai > A2. Since Ai, A2 are real we have 

which is only relevant if A + and A_ are complex conjugate of each other. A± are either 
both real for ( Al ~ A2 ) 2 ^ 5 2 or complex for 5 2 > ( Al ~ A2 ) 2 . In the last case we choose 



2 v v 2 

such that ImA + > 0. Now we want to integrate out and 5 for fixed A + ,A_ and find 
first 

<9(A+,A_) Ax-A 2 



9(Ai,A 2 ) A + -A_ 

which means 



(12) 



(Ai - X 2 )dX 1 dX 2 = (A+ - A_)dA + dA_ (13) 

valid in a sense of an alternating product of differentials, also in the case if A + , A_ are 
complex. Integrating out <fi and 5 we obtain for the measure of eigenvalues 

dfi(A+, A_) = — ^ r d</> / <W • c/A + rfA_ 2(A+ - A_) e -( A + +A - +4<52 )/ 2 (14) 
(2tt) 2 io Jrc 

where 5 is integrated over the region 1Z given by eq. (fTUj) . If A-t are real this yields 
simply 

d/x(A+, A_) = —^dA + dA_ (A + - A_)e-( A + +A -^ 2 (15) 
2y 27T 

and if A± are complex 

dfi(A+, A_) = -^dA + dA_ (A+ - A^)e- (A + +A - )/2 erfc(| ImA + |V2) (16) 
2v 27r 

with erfc(£) = ^= J z °° dxe~ x2 . Formula ffl6|) is in this form also valid for the real case. 
This measure is positive, since in the real case we assume A + > A_ and in the complex 
case A-t = x ± iy we have dA + cfA_ = (dx + idy)(dx — idy) = —2i dxdy and y > 0. 

We can also check the normalization. For the real eigenvalues we obtain 

dii{A+,A-) = ^=. (17) 
a± real V ^ 

This is the probability that both eigenvalues are real. For the complex part we obtain 

correctly 

rf/i(A + ,A_) = l-i= (18) 
A± complex V ^ 

such that total probability is 1. 
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2.2. General derivation 

We start again from eq. (CQ). The dimension N is even or odd, we will see that the odd 
case is more complicated. Again we will bring J by an orthogonal transformation O to 
a kind of lower triangular form 

J = O JO' 1 = 0(A + A) Cr 1 (19) 

with in the even case |iV 2x2 blocks A in the diagonal and in the odd case one lxl 
block together with |(AT — 1) 2x2 blocks in the diagonal. Then we obtain again 
eq. (j3J) for the increment dJ . Since 0~ l dO is skew-symmetric we can consider dJ and 
(0~ 1 dOJ — JO~ 1 dO) ahovc (i.e. where J is zero) as independent variables. The latter 
may be reduced to (0~ 1 dOh — AO~ 1 dO) abovc with Jacobian 1 due to the triangular 
structure of A, and these again to (0~ 1 dO) abovc with Jacobian rL<j(Ai — A?) wnere Aj 
are the eigenvalues of A. Here FJ means that the blocks in the diagonal are excluded: 
Yli<j = riabovo- If the blocks A are themselves diagonal this is easy to see. But it is then 
even generally true by diagonalization of the blocks A in the generic case. 

Thus we integrate first over O(N) /0{2) N / 2 in the even- A" case and over 
0(N)/0(2)( N ~ 1 ^ 2 x 0(1) in the odd- A" case, that is over the orthogonal group modulo 
rotations which leave the diagonal block structure invariant. Then we integrate over 
A. For the remaining blocks in the diagonal we can use our 2-dimensonal result from 
Section 12.11 Finally we arrive at the result derived by Lehmann and Sommers [3] and 
rediscovered by Edelman [I] 

d/i(Ai, A 2 , . . . , Atv) = K N ■ dAi . . . dA N ■ 

■ JI(A t - A,) ■ e-^/ 2 ■ (J] erfc(| ImA^)) 17 ' (20) 

i<j i 

with 

K N = VO{N) ■ 2- N {27t)- N{N+1)/A . (21) 

The constant K N comes from VO(N) /(4vr) Ar / 2 in the even case and VO(N) /(47r)( 7V ~ 1 )/ 2 -2 
in the odd case and integration over A. In both cases the result is the same. 

N 2ir d / 2 f\ 

^w=nf7dM=/n(°" ido )« ( 22 ) 

d=l V ' ' J i<j 

is the volume of the A^-dimensional orthogonal group. 

Here we have to assure that the transformation is unique. Thus the eigenvalues Aj 
have to be different and arranged in a special order. If all eigenvalues are real we 
assume A x > A 2 > A 3 .... If two eigenvalues are complex we assume Ai = A 2 and 
ImAi = — ImA 2 > and A3 > A4 > A5 . . .. If 4 eigenvalues are complex we assume 
Ai = A7, A 3 = Ai, ReAi > ReA 3 , ImAi, ImA 3 > and A 5 > A 6 . . . And so on. The 
measure is written in such a way that it is positive even in the complex case for which 
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Now we want to determine the correlation functions. To this end we go over to two- 
dimensional variables 

A — > z = x + iy with d 2 z = dxdy . (23) 

Then we ask for the symmetrized probability P(zi, Z2, ■ ■ ■ , Zf^)d 2 z\ . . . d 2 z^ that the 
complex eigenvalue tuple (z\, z 2 , ■ ■ ■ , z N ) appears in the volume element d 2 zi . . . d 2 ZN, 
where now we drop the restrictions for A, = Zi. 



3. Generating functional 

3.1. Definition 

We are interested in the n-point densities (or correlation functions) 

R n (z u z 2 ,..., z n ) = . . . j4-rZ\f\ (24) 

oj{zi) df(z n ) /sl 

with 

Z[f] = [ d 2 Z x . . . d 2 Z N P{ Zl , Zjj, ... , z N )f{ Zl ) . . . f{z N ) . (25) 



Recall that the variables Ak = Zk = Xk + iyk are considered here as two dimensional 
vectors. We obtain Z[f] by integrating over our joint density (|20l) . Here we have to 
distinguish between the different cases: all Zk real, 2 eigenvalues complex, 4 eigenvalues 
complex and so on: 

Z[f] =K N f d 2 Zl . . . d 2 z N f( Zl ) . . . f(z N ) Y[( Zi - Zj ) J] e"^/ 2 ( erfc(| lmz k \V2yf 2 

J i<j k 

■ {6(y 1 )6(y 2 ) ■ ■ ■ 5(y N ) 6(xi > x 2 > . . . > x N ) 

+(-2i)Q(y 1 )6(y 1 + y 2 )S(xi - x 2 )5(y 3 ) . . .5(y N )Q(x 3 > ...> x N ) 

+ (-2 i ) 2 &(y 1 )&(y 3 )5 2 (z 1 - z^)5 2 (z 3 - zi)G( Xl - x 3 ) 

■S(y5}S(y e ) ■ ■ ■ 5{y N )Q{x 5 > x 6 > . . . > x N ) + . . .} (26) 

By integrating over with the symmetric function f (z-\) . . . f (zjsf) the integrand will 
automatically be symmetrized. We used the notation 5 2 (z\ — z^) = 5(y\ +y 2 )S(xx —x 2 ). 
0(xi > x 2 > ■ ■ ■ > xn) stands for the restriction x\ > x 2 > . . . > xn- 

3.2. Vanderraonde determinant 

Let us write the Vandermonde determinant as 

llizi - zj) = {-ifi^V 2 Y[( Zl - Zj ) = (-1)^ det^-, z\~\ . . . , z^] (27) 

i<j i>j 
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with k = 1,2, ...,N. The determinant can be written with Grassmannian variables 



r)%,T]i and Berezin integration 



Y[(zi -Zj)= / drf x dr]i . . . dn* N dr] N exp (- ^ rf k z* \ ) ( 2> 

i>j kl 

Integrating out one set of Grassmannians (rji) we obtain: 

Y[( Zi - Zj ) = J dr)tdr,* 2 ...d V * N t^ 1 ) (]T t^" 1 ) . . . r^" 1 ) . (29) 

i<j ' k k k 

Note that the integrand factorizes in a product of identical functions of different 
arguments ( Zi ). 



3.3. Real case 

In the case that all eigenvalues are real we have to integrate a function f(x\) . . . f(xpf), 
where f(x) is actually Grassmannian from section [3T2l with the restriction 

9(xi - x 2 )Q(x 2 — x 3 ) ... <d(x N -2 - x N ^i)Q(x N ^! - x N ) (30) 

where Q(x) is the Heaviside step function = 1 for x > and for x < 0. 

Using Mehta's method of alternating variables we integrate first over x^, x^-2, %n-a 

and so on and obtain an integral 

■■■f(xN-3) / dx N - 2 f(xN-2)f(xN-l) / dx N f(x N ). (31) 

Jxn—i J —OO 

Now we use that f(x) is Grassmannian and therefore also dx^fix^) for which the 
square vanishes. Thus we may replace the above result by 

...f(x N - 3 ) / dx N ^ 2 f{xN-2)f{xN-i) I dx N f(x N ). (32) 

J —oo J —oo 

There remains the restriction xi > x 3 > . . . > xn-i for even N, dropping it we obtain 



(N/2) 

And if N is odd we have 

"+00 

1=1 dxf(x 



+°° ~ /•si „ -| AT/2 

da:i/(xi) / dx 2 f(x 2 ) 



((N-l)/2)\ 



+00 



dxif(xi) / dx 2 f(x 2 ) 



1 (JV-l)/2 



(33) 



.(34) 



5.^. General case 

Now we first assume that iV is even and all eigenvalues are pairwise complex conjugate. 
We have the restriction (besides Z\ = z^, z 3 = etc; yi, y 3 , ... , yjv-i > 0) 

Q(xt - x 3 )Q(x 5 - x 7 ) . . . 6(xat_ 3 - x N -i) . (35) 

Dropping the restriction we obtain for the integral over f(zi) . . . f(zm) 



(AT/2)! 



2i / d 2 zQ{y)f{z)f{z) 



N/2 



(36) 
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Now it is easy to sum over all mixed cases. For even N we obtain 



I = 



(JV ^|J(T)[- 2i / d2 ^w^ 



M 



dxif(x l ) I dx 2 f(x 2 ) 

oo J — oo 



(N/2)! 



-2i J d 2 ze(y)f( 



z)f(z) + 



"+°° ~ r x i ~ -i n/2 
I dx 1 f(x 1 ) l dx 2 f(x 2 ) 



and similarly in the odd case we obtain 

"+°° „_ i 
dxfix)-^- 

> V 2 >■ 



-2i / rf 2 z9(|/)/(z)/(z) + 



/+°° ~ r x i ~ -I (jv-i)/2 
dxif(xi) I dx 2 f(x 2 ) 
oo J —oo 



(37) 



(38) 



5.5. Z[/] as a Pfaffian 

Now we may write the generating functional as integral over Grassmannians (in the 
even-iV case) 

AT/2 



A? 



c??72 . . . d?7 



iV 



(iV/2)! 



1..JV 



kl 



= K N W.. rf *exp(-lE^*) 



with the skew-symmetric matrix 

4u = y rf^ 1 rf^2/(^l)/(^2)^(^1^2)^- 1 4" 1 

and the skew-symmetric measure 

T{z x , z 2 ) = e -rt + ®' 2 \2i5\ Zl - z 2 ){Q( yi ) erfc(y 1 v / 2) 



(39) 



(40) 



- 9(y 2 ) eric(y 2 V2)} + 5{ yi )5{y 2 ) (Q(x 2 - x ± ) - G( Xl - x 2 )) 



(41) 



We have antisymmetrized because rf k are Grassmannians. Such an expression for Z[f] 
is called a Pfaffian 



Z[f] = K N Pfaff(A) = AWdet(A) . 



(42) 



The Pfaffian is an analytic square root of the determinant det(A). It is only defined 
for an antisymmetric matrix. In this case we know that for positive f(z) P faff (A) is 
positive, such that the square root of the positive determinant det(A) is also positive. 
We actually need f(z) only in an infinitesimal region near f(z) = 1. 
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We may immediately consider odd N. In that case Am has always a zero eigenvalue 
and cannot be inverted (what we will need for perturbation expansions). Therefore we 
increase artificially the number of Grassmannians by 1 



Z[f} = K N / dr)l...dr]* N+1 r]* N+1 / dxf{x)e '"' 2 ^ '/ v 



N 



k=l 



nm 

Now n, m = 1, 2, . . . , JV + 1 and 



nm'lm 



Am Ck 
-Ci 



(43) 



(44) 



with 



C fc = y dxf(x)e~ x2/2 x k - 1 . (45) 

Thus we obtain in the odd-iV case 

Z[f] = K N Pfaff(5 nm ) (46) 

The correlation functions can be found by multiple derivating of Z[f] w.r.t. f(z) at 
/(*) = I- 

4. 1-point density 



To calculate the 1-point density we use 



Rx(z) 



5)nZ[f] 



5f{z) 



(47) 



f=i 



4-1. Even N , complex eigenvalues 
In the even-iV case we have 



ltflndetA _ 1 1 5 A 
l[Z) " 2^JW~ ~ % A Jf(z) 



(4* 



where Am is given by eq. ( HD1 . Let us call Afcj = Am\j=i and introduce the kernel 

1...N 

JC N (z 2 ,Z 1 ) = J2 A kl Z 2~ lz t 1 - ( 49 ) 
fc,/ 

Then we obtain 

(50) 
(51) 



Ri{zi) = j d 2 z 2 J r (z 1 ,z 2 )JC N (z 2 ,z l ) 

with T{z\, z 2 ) given by eq. (j4ip . If we insert eq. (j4~Tl) we obtain 
i2i(zi)=i2f(«i)+«J(yi)i2f(x 1 ) 
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B°{z x ) = 2i sgn ( Vl ) erfc(\y 1 \V2)e- x2 i + rt)C N (zI,z 1 ) 

and 

r+oo 

R?{xi) = / dx 2 sgn (x 2 - Xi)e" {a; ? +x 2 )/2 /C Ar (x2,Xi) 
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(52) 



(53) 



which is a smooth part Rf (z\) in the complex plane and a part i?f (x\) concentrated 
on the real axis. 



If we compare i?f (z\) with Edelman's expression [lj for the complex 1-point density we 
find for JCn(z 2 , Z\) using that z~i and z\ are independent variables 



JCjv(z 2 , z l) 



AT-2 , , 
z 2 - Zisr^ ( z l z 2) 



1...N 



E 



2V27T rai 

n=0 



(54) 



Thus surprisingly, the skew symmetric matrix A fcZ has a very simple tridiagonal 
structure 

" -I 



A ki 



1 



2V27T 















1 



(JV-2)! 



(AT-2)! 




(55) 



which leads to Edelman's expression 

2|y 



AT-2 



erfc (|y|V2)e- a;2+2 ' 2 £ 



\2n 



n=0 



7V. 



(56) 



For N — 2 it agrees with our 2-dimensional expression eq. f|T6|) . Note that i?i(2;) as a 
density is normalized to N: J d 2 zR x (z) = N. Similarly / d 2 Z\d 2 z 2 R 2 (z\, z 2 ) = N(N — 1) 
and so on. 
Using the formula 



N 

E 

n=0 



n! 



JV 



e" u 



u 

ivT 



we can also write 



2|y 



AT-2 



erfc (V^lyDe* / ^e-- 

27T J\ z \2 1 (iV - 1, 



(57) 



(5f 



analytic in N. From fl55l) it is easy to check the normalization (eqs. (121I22P ) using the 
duplication formula for the Gamma function. Note that Pfaff (Am) is positive. 
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4-2. Even N , real eigenvalues 

Since the same analytic kernel appears also in the density of real eigenvalues eq. (l53~I) 
we have using eq. ( 1571) 

1 



2V27T 

Integration by parts yields 



" ■ i ■ > ■- / °° u N ~ 2 

dx 2 sgn (x 2 — Xi)e 5 [x 2 — xi) I due~ u - 



(JV-2)! 



(59) 



^y |!Bl p (iV-2)! 
j r+oc 



2V27T 7_ 



dx 2 sgn (x 2 - a; 1 )^^C^e-^D/ 2 



(iV-2)! 



with 



Rf( Xl ) - D N x^-\- x ' /2 Q{N odd) 



I due U -— T — — + 

+ v^X r(iv-i) 



(60) 



|x 1 | JV - 1 e- ai ?/ 2 



and 



iV 



dxe x2 l 2 x N 2 



(61) 



(62) 



, V2^T(N - 1) 2 N / 2 T(N/2) ' 

The second term in the last line of eq. (|60|) appears formally only if iV is odd. We 
see that -Rf(xi) is an analytic function of N which is manifestly positive. This result 
has been obtained by Edelman, Kostlan and Shub [5] and is also valid for odd N. We 
have to subtract just the second term in flBDl) which is proportional to 0(iV odd) (which 
means that this is only 1 if N is odd otherwise it is zero) to get the correct answer valid 
for even and odd N. We will see the consequences in the following. 



4.3. OddN 

If N is odd R\{zx) is given by 



Ri(zi) 



1 Tr 1 SB 

2 B6f{ Zl 



1...N 



k,l 



N 



&n+i,i 



i=i 



f=i 



f N 

/ d 2 z 2 J r (z 1 , z 2 )!C N (z 2 , Zl ) + B N 1 +ij s (Vi) e ~ 3 * /2x i' 1 ■ 
J i=i 



(63) 
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Since the complex-eigenvalue part Rf(zi) is again given by Edelman's expression (|5"6"1) 
and the real-eigenvalue part i?f (xi) is given by eq. flBTj) we conclude that B~^ has again 
a very simple tridiagonal structure 





B 



-i 



2V27T 





0! 



0! 







o h 







1 



{N-2)\ 





(AT-2)! 

2^2ttD n 





-2y/27rD N 




(64) 



only the last row and column make the matrix B 1 invertible. The kernel now is given 
by 



1..JV 

K N (z 2 ,z 1 ) = J2B k }zt 1 z[- 1 

k,l 



Z 2 ~ Z X 



N-2 



\ Z 2 Z \ j 



2V27T *-^L n\ 

n=0 



while 



D— 1 __R-1 _A - 
n N+l,k - n k,N+l - °k,N 2 iV/2r(jV/2) 



N 



(65) 



(66) 



In a diagrammatic representation, which we will use later, we have with an obvious 
notation 1 1 



Ri(zi) 





x N-l e -4/2 



^(1,2)^(2,1) + %)^^ 



6(A^ odd). 



(67) 



Note that the last term gives exactly the correct density for iV = 1. 



5. Correlation functions 

5.1. Diagram expansion 

Let us write f(z) = 1 + u(z) and 

A kl = A kl + C kl = J d 2 Zl d 2 z 2 z\- x z l 2 - x T{z x , z 2 )(l + u(zi))(l + u{z 2 )) (68) 

and try to expand Z[l + u] in powers of u(z) to get the n-point densities. To this end 
we first expand Z[l + u] in powers of C k \ (first the even-iV case): 

Z [1+Ml ^„/^...^e Xp (-I D . H + Cw) «) 

J kl 
n=0 ' \ kl 



(69) 
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We use now that there is a fermionic Wick theorem with 

(vM)o = All = k • • I (70) 

-C kl = (71) 
such that with the linked cluster theorem 

In(Z[l+«]/Z )= (*\ { \ { \ (72) 



or 



Z[l+u}/Z = 1 + 




(73) 



which includes for each internal point a summation, a minus sign for each closed fermion 
loop and the factor 1/ (order of invariance group) for each diagram. It is important for 
the symmetry, that in the expansion of the Pfaffian the fermion lines carry no direction. 
However, translating the diagrams one has to go through it in a certain direction, which 
determines the sign of each element. If we went through it in the opposite direction, all 
elements would obtain the opposite sign, because they are skew-symmetric. The result 
is the same. Z Q is just given by Z = Z[l] = 1. 

Introducing the kernel K.n(zi,z 2 ) from eq. ( 1541) and T{zi,z 2 ) from eq. ( 1411 we 
reinterpret the diagrams as 

JC N { Zl ,z 2 )=1C N {l,2) = l* .2 (74) 

- T{z\, z 2 )(u(z 1 ) + u(z 2 ) + u( Zl )u(z 2 )) = 1 • • 2 (75) 

and at each internal point z we have now an integration over d 2 z. In an obvious notation 
we now produce all correlation functions by functional derivatives w.r.t. u(z): 

^6+{do + dwdoo ■■■}+••• < 76) 



1 2 



5 2 Z[l + u] 
Su(zi)5u(z 2 ) 



OIQOQQ 



+ 



(77) 



OQO 



+ 



(7f 



General Eigenvalue Correlations for the Real Ginibre Ensemble 



14 



And so on. A label 1 means a functional derivative w.r.t. u(zi) and thus at this vertex 
appears no integration. At a wavy line one can only differentiate twice and only at 
different vertices, otherwise there is no contribution. To obtain all correlation functions 
or n-point densities we have to put u(z) = at the end. Then a lot of diagrams 
disappear. In the following diagrams we reinterpret a wavy line as: 

l*~~~v. 2 ^-J r (z 1 ,z 2 ) (79) 

and obtain 

1 



Ri(zi) = 




«,<*,*) =q + 00 + O + O 2 (8i) 

2 2 
And so on. These diagrams have no longer an invariance group, because the external 
vertices, which are not integrated over are labeled and distinguishable. Here no wavy 
line is possible that has only internal vertices. Differentiating a diagram, which contains 
( 1751) . once w.r.t. u{z\) at u = yields a wavy line with one external vertex (1). Differ- 
entiating a diagram, which contains (!75j) . w.r.t. u{z\) and u(z 2 ) at u = yields a wavy 
line with two external vertices (1,2). 

We could also calculate the cluster functions which are in each order only the connected 
diagrams, i.e. here the one-loop diagrams. 

We will see that in the case of odd iV the diagrams are slightly modified. There is an 
additional graphical element for the correlation functions (independent of the direction 
going through) 

7 N ~ l 

1 2 = Hpjjpj • e-^5(y 2 )Q(N odd) (82) 

The cross corresponds to the additional artificial Grassmannian and therefore can only 
appear once in each diagram. The result for the correlation functions is then 

1 1 

/?,(:.,.) = {} + () m) 



R2(z 1 ,z 2 ) 



0*0 

6*66*66*66 
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+ 



0*0*0 

2 2 2 

2 

+i O +i O + O 1 (84) 



And so on. The terms with the crosses, which are only present if N is odd, are obtained 
by labeling in all possible ways one internal vertex in the original diagrams by a cross. 
One has still to keep in mind that closing a loop yields an additional minus sign. 

5.2. Correlation functions, N even 

Let us translate the diagrams again (for even N): eq. (l80l) leads to eq. ( |50l) and eq. 
(JHU) to 



R 2 (zi,z 2 ) = J r (z 1 ,z 2 )IC N (z 2 ,z 1 ) + j d 2 z 3 d 2 z 4 { 

+ J r (z 1 , z 3 ) K, N (z 3 , z x ) T{z 2 , Z4) /Ctv0 4 , z 2 ) 

- T(z x , z 3 ) K, N (z 3 , z 2 ) T(z 2 , z 4 ) JC N (z 4 , z x ) 

- F(zx, z 3 ) fC N (z 3 , z 4 ) T{z^ z 2 ) JC N (z 2 , zi)} (85) 

And so on. Now we know how to generate general n-point densities. We observe that at 
the external vertices there are no integrations. Thus the diagrams are cut into factors 
which are special diagrammatic elements. These are used to build for an n-point density 
closed loops with n external vertices, using the rule that a wavy line can only be linked 
with a straight line, as is done in eq. (1841) . In the following we will show that the 
result is just what is called a quaternion determinant of a self-dual 2n x In matrix 
(M = l,2,...,n): 



Rn(zi, 22, ... , z n ) = (-1)" qdet 



( 1 A 

k* • »l k* •/ 



(k • • / + k I ) k I j 



(86) 



The entries are the appearing graphical elements. The (— l) n reminds us that we have to 
take into account a minus sign if we close a loop. The sign of the quaternion determinant 
is defined in such a way that in the expansion of R n (zi, z 2 , . . . , z n ) there appears the 
positive term R\(zi) ■ Ri(z 2 ) ■ . . . ■ Ri(z n ) coming from the diagonal elements, which 
gives the behaviour for large separation. We will see that the quaternion determinant 
is related to a Pfaffian. Since by construction integration of R n over d 2 z n leads to 
(N — n + l)R n _i, eq. fl86|) implies an integration theorem [6] for this type of Pfaffians. 
This however includes here all 5-type contributions to R n (-5-functions for real eigenvalues 
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and complex conjugate pairs). To separate all these terms may still require some 
combinatorial analysis. 

5. 3. Expansion of a Pfaffian 
Let us consider 

/ 1...2M \ 

Z{e) = J d V t... d V * 2M exp f -- J2 vKAh + CM J (87) 

now with Aki = eJki and J = ) with M x M entries. Obviously Pfaff (Cm) = 
lim £ ^ Z(e). Now we expand Z(e) in powers of C using the fermionic Wick theorem with 
A k l = —\Jkh Since Z = (-l) M ( M - 1 )/ 2 £ M in the limit e — > only terms with power M 
of A -1 survive, higher powers do not occur. On the other hand (^77* )o = — \Jh- Thus 
we have 

(_-qm(m-i)/2 pfaff ^ = gum of all diagrams of order M 




o 00 o 



+ . . . 

M times (M - 2) times 
= (-1) M qdet (J • C) = (-1) M qdet ( • • • ) (88) 

with • • = — J,»v~~v~v^« — —C. 

We see that only the matrix JC appears which is selfdual (i.e. J(JC) T J T = JC 
if C = — C T ). We also may replace the expansion of the Pfaffian by an indexed 
diagram expansion in which all external vertices are different. To see this formally 
let C H -> X k C k iXi with X k = 1 for k = M + 1,M + 2,...,2M and differentiate 
Pfaff (ACA) w.r.t. Ai, A2, • • • , Am- This produces the indexed diagram expansion and 
since Pfaff (ACA) = Pfaff (C) • A1A2 • • • \m the result is the same. This proves our claim 
( I86|) . Note that if one uses a different definition of the Pfaffian with a different order of 
Grassmannians the factor (— l) M ( A/ ~ 1 )/ 2 in (1881) may be canceled. 

5.4- Correlation functions, N odd 

Let us recall for odd N the diagrammatic expansion on the level where the vertices carry 
the number k of Grassmannian n* k . Then again 

(«)o = ^ 1= A:« •/ 

-C kl = k • • I (89) 

with k, I — 1, 2, . . . , N and there are additional graphical elements 
(v*N+iV*k)o = B N \ lk = x • k 
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alxe- x2/2 x l u(x) = I* x (90) 



We have from eq. flMJ B^\ l k = 5 kN /2 N / 2 T{N/2) = -B^ N+1 and find that in the 
diagrammatic expansion only the combination 

k ^—l =- 2 J^ N/2 y J dxe-* 2 ' 2 x l - l u{x) (91) 

occurs. In this diagram it is not important in which direction one goes through it: the 
result is the same. Thus one obtains the additional diagrams: 

00 

Only diagrams with one cross are needed. At each vertex without cross there appears 
a /c-summation over z k ~ 1 . 




o 




. . . (92) 



Now we again may reinterpret the diagrams as: 
!• »2 =/0v(l,2) 

1 . . 2 = - 2)(u(l) + u(2) + u(l)u(2)) 

1. — * — .2 = -u(l)e- x * /2 5(y 1 )z?- 1 /2 N/2 T(N/2) (93) 

with integration at internal vertices. Remember that a closed loop yields a factor (—1). 
Finally we differentiate the diagrams w.r.t. u(z) several times and put then u(z) = to 
obtain all correlation functions. Note that there is no differentiation possible at a cross. 
The result is the correlations fl83l |84|) and corresponding higher orders. 



The general formula for the n-point densities may again be written as a quaternion 



determinant with k, I = 1, 2, . . . , n: 



R n (z u . . . , z n ) = (-l) n qdet 



' (k i + k •— • /) k* •/ 

^ k mssssssml [k I + k + ™K • I) 



(94) 



where a label k means the site Zk and 

k mxssw I = k • »l + k I + k X— »~« I + k I (95) 

Here the diagrams are interpreted as in eqs. (1741) . f[79|) . (1521) . Since the quaternion 
determinant is a Pfaffian and the crossed terms factorize, it is easy to see due to r] 2 = 
for a Grassmannian 77, that only the first order terms in an expansion of powers of the 
cross contribute. We also have not to worry about the 5-functions which appear inside 

the quaternion determinant in all elements except • • , since we see from the 

expansion of the Pfaffian that each 5-function appears at most once at each site. The 
5-functions single out special correlations of real eigenvalues (8{y^)) or pairs of complex 
conjugate eigenvalues (5 2 (zi — zj)). 
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R-n\Z\i • • • i %n) 



-l) n qdet 



/ G(l,l) ■ 


■ G(l,n) 


K(l,l) ■ 


•• K(l, 


n) \ 


G(n,l) • 


■ G(n,n) 


K(n,l) ■ 


•• K(n. 


n) 


W(l,l) ■ 


■ W(l,n) 


G(l,l) • 


■■ G(n. 


1) 


\W(n,l) ■ 


■ W(n,n) 


G(l,n) • 


■■ G(n, 


n) / 



(96) 



with 

G(l, 2) = 1 • — • — • 2 + 1 • — * — • 2 

K(l, 2) = K N (1, 2) = 1 • • 2 

W(l,2) = l« »2+l« • — • — •2 + 1 • — * — • — «2 + l« — • — * — »2 (97) 

The quaternion determinant of a matrix C is according to eq. (188]) up to an overall sign 
equal to the Pfaffian of J ■ C. The sign is in our case easy to find since all R n ^ 0. 
G(l, 2) may be split into two parts depending on wether the second argument is real or 
complex: 

G(1,2) = G C (1,2)+G*(1,2)<% 2 ) (98) 
Then we obtain 

G c (l,2) =/Cjv(zi,^)2isgn {y 2 )e 



-*l+vl er fc(|y 2 |V2) 



u 



N-2 



due , 
2v2tt J Z1 .^ (N-2)\ 

Hi 



2i sgn (y 2 )e 



-l+vl er fc(|2/ 2 |v^) (99) 



analytic in N. In the second part G (1, 2) we can make a partial integration like in 
section 14.21 which cancels the term with the cross, which is present only for odd N, and 
obtain 



G R (1,2) 



1 



2ll 
1 



u 



N-2 



due u - 

(N - 2) 



X-2 



dxe x2/2 



N-2 



JV-1 



'2vr Jo (iV-2)! 1 

which is analytic in N and valid for even and odd N. 

Finally we consider W(l, 2) which splits into 4 parts 

W{1, 2) = W cc (l, 2) + <J(y 1 )W flC (l, 2) + 11^(1, 2)<% 2 ) 
+ <% 1 )% 2 )W iiR (l,2) 



(100) 



(101) 
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First we obtain W cc (l,2): 

W cc (l,2) = -2i8 S n(y 1 )e-^eiic(\y 1 \y/2){g 2 (z 1 -z5) 
+ TC N {zi,zl)2i sgn(y 2 )e-^+% 2 e rfc(|y 2 | v^)} 

= - 2i sgn( yi )e- x ' + y' erfc(| yi | V2) {S 2 ( Zl + G c (zT, z 2 )} . (102) 

Again in W RC (1, 2) and W CR (1, 2) we can make a partial integration to cancel the cross 
term and obtain 

W RC (1,2) = -W CR {2, 1) = 2i sgn{y 2 )e-^ +y " erfc(|y 2 |v^)G fl (^,a;i) (103) 
again valid for even and odd iV and analytic in N. 

Finally the most complicated term W RR (1,2) can also be reduced to G R (1,2). One 
splits off a factor e~( x l +x l)/ 2 and derives a first order differential equation in x±. Using 
the skew-symmetry and again a partial integration to cancel the cross-terms one arrives 
at 

W RR (1, 2) = 2 f 1 dx e~^ +x2 ^ 2 G R (x, x 2 ) - e ~^ +x ^ 2 sgn(x 2 - x ± ) 
= -/ dx e-^ +x2 )/ 2 G fl (a;,a;i)+e-( x ? +x2 )/ 2 G fl (a;,a; 2 ) 

J Xl 

- e -^+^)/ 2 sgn (a; 2 - Xl ) . (104) 



Again this expression is valid for even and odd N and analytic in N. The second line 
is a nontrivial consequence of the skew-symmetry of ^^(1,2). In all cases the ana- 
lyticity in N is easily seen for arguments zi, z 2 , . . . , z n positive, but then can be extended. 

At the end of this section let us write down some correlations in the notations of this 
section: 

=-G(l,l) (105) 
R 2 (l, 2) = G(l, 1)G(2, 2) - G(l, 2)G(2, 1) - W(l, 2)K(2, 1) (106) 

The last two terms yield the connected part, i.e. the cluster function 

R c 2 on {l, 2) = -G(l, 2)G(2, 1) - W(l, 2)K(2, 1) . (107) 

In general one draws all possible diagrams with elements G, K, W and the sign 
^_ ^number of fermion loops 

5.6. Numerical evaluation 

Let us introduce the function <j)(z): 

<j>{z) = -2% sgn(y) e~ x2+y2 erfc(|y| y/2) = -</>(z) (108) 
and the incomplete Gamma function 7*(n, x): 

T(n)x n 1 *(n,x)= f due^u 71 ' 1 . (109) 
Jo 
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7* (n, x) is an analytic function of both arguments. It has the power expansion 

oo m 

7 *(n,x) =e- x V — -. (110) 

^ T(n + m + l) 



m=0 



Then we obtain as basic functions from eqs. (|54l) . (!571) . ( 1991) . f llOOp 



if (1,2) = K N {z u z 2 ) = Z ±-^^{\ - (z 1 z 2 ) N - 1 1 *(N-l,z 1 z 2 )) (111) 

2v 2tt 



and 



and 



G c (l,2) = G c ( Zl ,z 2 ) = -Kxiz^zq)^) (112) 



G R (l,2) = G R (z u x 2 ) = — -Le-^+^l - {zix^-^^N - l, Zl x 2 )) 

V 2n 

e- x y 2 ( Zl x 2 ) N - 1 ,/iV-l xl 



7*( 1 n ) ■ (H3) 



2*-V2r(JV/2) 
From this we obtain 

W cc (l,2) = 0(^)(«5 a (zi -2*) + G c (-zT,z 2 )) 

= <t>{z 1 )5\z 1 -zi) + W cc {1,2) (114) 

The first term leads to a self-correlation of a complex conjugate pair. Then 

W RC \l, 2) = -jy CK (2, 1) = -<P{z 2 )G R {- 2 , x x ) (115) 

and W RR (1,2) is given by eq. (11041) . These formulae immediately yield all n-point 
densities: 

R c ; (1) = -G c (l, 1) = K N {z lt zO<l>{zi) (116) 
J2?(l) = -G R (1, 1) = -1,(1 - ^"V(JV - 1, x?)) + 

V Z7T 



e ^i/ 2 ^^ ^ *(N — 1 x 2 \ 

+ 2 jv -v2r(A^/2) 7 I 2 ' ^ ^ 

and also ^(1,2) splits into 5 parts: 

i? 2 (l, 2) = Rf (l)5\ Zl -zi) + R% c (l, 2) + R RC (1, 2)8( yi ) + R^ R (1, 2)8(y 2 ) 

+ R RR (l,2)8( yi )S(y 2 ) . (118) 

All the terms follow from eq. (I106p . The first term corresponds to the correlation of a 
complex eigenvalue with its complex conjugate. The other smooth terms correspond to 
correlation complex-complex, real-complex, complex-real, real-real. For completeness 
let us write down these terms: 



R 





:i,2) 


= G°\ 


;i,l)G c (2,2) 


-G c {l,2)G c {2,\) 


- w cc 


(1,2)^(2,1) 




• RC [ 


;i,2) 


= G R [ 


1,1)G C (2,2) 


-G c (l,2)G fl (2,l) 


- W RC { 


;i,2)if(2,i) 






;i,2) 


= G C [ 


1,1)G*(2,2) 


-G c (2, 1)G R (1, 2) 


-W CR ( 


;i,2)if(2,i) 






:i,2) 


= G R [ 


1,1)G R (2,2) 


-G*(l,2)G fl (2,l) 


- W RR ( 


;i,2)/f(2,i) 


(119) 
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Figure 1. 1-point density R%(zi) from simulation for N — 20. The total integral is 
normalized to N. 




Figure 2. Left: 2-point density i22(2i, 2a) as function of z\, for fixed = 2 + 2i and 
the same simulation data as in Fig. [T]with N — 20. The total integral is normalized 
to (N — 1) • i?i(z 2 ). Right: The same for a simulation with N = 15 and z 2 = 2 + 0.5i. 



In the following we do some numerical simulations and compare them with numerical 
evaluations of the above formulae. We draw randomly matrices from the Gaussian 
ensemble ([[]) and plot at first a histogram for the eigenvalues in the complex plane, 
which yields R\{z\) (Fig. CD). We see that the eigenvalues lie in a circle with radius of 
order \fN and that a finite fraction lies strictly on the real axis repelling the remaining 
pairs of complex conjugate eigenvalues from the real axis. Then we take the same set 
of eigenvalues and choose only that subset with one eigenvalue close to a fixed value Z2- 
Plotting a histogram of this set we obtain ^2(^1, z 2 ) (Fig. [2]). We see that the complex 
eigenvalue z% repells all the other with a cubic law in distance and that there is again a 
finite fraction of eigenvalues on the real axis corresponding to complex-real correlation. 
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Figure 3. 3-point density Ra(zi, Z2, Z3) from simulation for N 
Z2 = 2 + 2i, Z3 = —2 + 2i as function of z\ (unnormalized) . 



20 and fixed 



0.1 



Figure 4. Analytical result for R^(xi) for N = 20 as a function of x% from eq. (JTT7J) . 



Furthermore one finds a 5-peak at the complex conjugate site I2 and also repulsion from 
that point. We can even take the same data, fix two eigenvalues z 2 and z 3 and plot 
R 3 (zi, z 2 , z 3 ) as function of Z\ (Fig. [3]). 

Using the above formulae we can calculate i?f(l), #f(l), R% c (l,2), R* c (l,2), 
R2 R (l,2) etc. exactly. Below we plot the functions i?f(xi), R^(zi), R^°(xi, z 2 ) and 
R§ c (zi, Z2) for the same fixed z 2 as in the simulation (Fig. HE). A further plot (Fig. 
|8]) shows the correlation R2 R (x, 0) for N = 2,3, 4, 5, 6, for which we have the analytical 
formula from eqs. ffTT9l [Till fTM flTTD : 



R RR (x,o) 



' R R (x)-—e 
2tt 11 ^ 2vr 



+ 



x 



-x 2 /2 



X 



erfc(^P 



(120) 



2V27T "^2 

The only iV-dependence is in the first term R R (x), which is almost constant in x 
for sufficiently large N. For small \x\ we have linear level repulsion. The connected 

2 

correlation at the special point is independent of N and decays with e~ x . 
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Figure 5. Analytical result for Rf(zi) for N — 20 as a function of z\ from eq. (|116p . 




-7.5 -5 -2.5 2.5 5 7.5 -7.5 -5 -2.5 2.5 5 7.5 



Figure 6. Left: Analytical result for the correlation R^ c (xi, 22) for N = 20 and fixed 
Z2 = 2 + 2i as a function of the real x\ from eq. (|119[) . Right: The same for N = 15 
and z 2 = 2 + O.bi. 

6. Conclusions 

In this paper we have derived, starting from the Gaussian ensemble, closed analytical 
expressions for all correlation functions, i.e. n-point densities, of eigenvalues of the real 
Ginibre ensemble of real asymmetric matrices, which is invariant under real orthogonal 
transformations. These n-point densities are not absolutely continuous. They contain 
contributions which are concentrated by 5-functions on the real axis and also point- 
measures for pairs of complex conjugate eigenvalues which are always present, because 
the eigenvalues of a real asymmetric matrix are either real or pairwise complex conju- 
gate. The n-point densities are written as quaternion determinants of certain 2n x 2n 
selfdual matrices, or as Pfaffians of corresponding skew-symmetric matrices. The Pfaffi- 
ans can be derived from a zero- dimensional fermion field theory, similar in structure as 
the matrix Green functions in the Nambu space for superconductivity All the n-point 
densities are expressible in terms of a skew-symmetric measure F{z-y, Z2) containing one 
part concentrated on the complex conjugate pairs Z\ = Z2 and one part on the real 
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Figure 7. Left: Analytical result for the correlation R 2 C (zi, z 2 ) for N = 20 and fixed 
z-2 = 2 + 2i as a function of z\ from eq. (|119p . Right: The same for N = 15 and 
z 2 = 2 + 0.5i. 




-6-4-2 2 4 6 



Figure 8. Analytical result for the correlation R^(x, 0) for TV = 2, 3, 4, 5, 6 from eq. 

poll . 

eigenvalues z\ = z~i, z 2 = ~z~2. For the N- dimensional matrix one constructs from 
T{z\,Z2) a skew-symmetric kernel £#(21,22) which together with .F^i,^) yields the 
building blocks for the correlations. 

To calculate £#(21,22) one has to invert an iV x iV antisymmetric matrix Ajy related 
to T[z\,z<i), which looks complicated, which however turns out to yield a very simple 
tridiagonal structure for Aj~j . To find this it is enough to compare the Edelman result 
for the complex 1-point density with the general form of the 1-point density. This is 
sufficient to obtain all correlations in the case of even dimension N. In the case of odd 
dimension N one has to increase the dimension artificially by 1 and has to invert instead 
the (N + 1) x (N + 1) dimensional skew-symmetric matrix B nm . Again a simple argu- 
ment using the result of Edelman, Kostlan and Shub for the real 1-point density reveals 
that has again a simple tridiagonal (but slightly more complicated) structure. In 
this paper we do not discuss detailed asymptotics for large N, which has partly been 
discussed elsewhere [H EJ |9] and which follows essentially from the asymptotic kernel 
£jv(2x, Z2) — (zi — Z2)e +ZlZ2 j2\phx. However we show that the different formulae of the 
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correlation functions for even and odd N can be combined always in one formula which 
depends on continuous and even analytic functions of N. 

Finally we have presented some numerical simulations which make clear the complicated 
structure of the correlations, for example the 2-point correlation. There is a smooth 
background of complex-complex correlations, then a point measure for a complex con- 
jugate pair, a part with one 5-function concentrated on the real axis for complex-real 
correlations and a part corresponding to real-real correlations doubly concentrated on 
the real axis. In comparison we have shown also some numerical evaluations of the 
analytical formulae. 

For the future there remain to be discussed more detailed asymptotics and more 
sophisticated functions like level spacings and distributions of extremal eigenvalues. 
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